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Abstract 

Optical and magneto-optical properties of YVO3 single crystal were studied in FIR, visible, and 
UV regions. Two structural phase transitions at 75 K and 200 K were observed and established to 
be of the first and second order, respectively. The lattice has an orthorhombic Pbnm symmetry 
both above 200 K as well as below 75 K, and is found to be dimerized monoclinic Pbll in between. 
We identify YVO3 as a Mott-Hubbard insulator with the optical gap of 1.6 eV. The electronic 
excitations in the visible spectrum are determined by three d-bands at 1.8, 2.4, and 3.3 eV, followed 
by the charge-transfer transitions at about 4 eV. The observed structure is in good agreement with 
LSDA+{7 band structure calculations. By using ligand field considerations, we assigned these bands 
to the transitions to the ^^2g, ^Eg Tig, and ^T2g states. The strong temperature dependence of 
these bands is in agreement with the formation of orbital order. Despite the small net magnetic 
moment of 0.01 fiB per vanadium, the Kerr effect of the order of 0.01° was observed for all three 
d-bands in the magnetically ordered phase Tn^cI < 116-?^. A surprisingly strong enhancement of 
the Kerr effect was found below 75 K, reaching a maximum of 0.1°. The effect is ascribed to the 
non- vanishing net orbital magnetic moment. 

PACS numbers: 71.20.-b, 78.20.-e, 78.30.-j 
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I. INTRODUCTION 



The interplay between kinetic energy and Coulomb interaction in strongly correlated 
electron systems enables small external or internal perturbations to manipulate electronic 
and magnetic properties of the matter. If, in addition, there is a ground state degeneracy, 
further interesting phenomena, like an orbital ordering can arise^. 

In yttrium orthovanadate, YVO3, spin and orbital ordering was suggested to be respon- 
sible for at least one of the temperature induced sign reversals of magnetization-. In the 
magnetically ordered phase, YVO3 is a canted antiferromagnet, where a net ferromagnetic 
moment is formed by a small angle between the antiferromagnetically (AFM) oriented spins 
of ions. The transition from the room temperature paramagnetic phase to the antifer- 
romagnetic state occurs at the Neel temperature T/v = 116 K. Upon further lowering the 
temperature, the magnetic moment gradually changes sign around 90 K. It has been argued 
that this magnetization reversal is due to the opposing effects of the Dzyaloshinsky-Moriya 
interaction and the single-ion magnetic anisotropy on the spin canting directioit^. The mag- 
netization switches sign again around 77 K. At this temperature the antiferromagnetic order 
changes from the high temperature C type to the low temperature G type-. A C-type AFM 
order corresponds to an AFM arrangement in planes and a ferromagnetic arrangement be- 
tween the planes. A G-type AFM order implies completely antiferromagnetic arrangement, 
both within planes and between them. Ren et. al?-^^ argued that along with the spin order 
there is an orbital order on the vanadium sites. According to the proposed picture, one of 
two 3(i electrons always occupies the xy orbital, while another electron occupies alterna- 
tively either the xz or yz orbital on different atoms, thus forming orbital order. In order to 
minimize the exchange energy a C-type spin order is accompanied by a G-type orbital order 
and vise versa. The onset temperature of the orbital order was reported to be much higher 
than the Neel temperature, Tqo = 200 K— . 

The proposed picture of the electronic structure still requires confirmation from, for in- 
stance, spectroscopic experiments. The strength of optical transitions involving a charge 
transfer between orbitals of different ions should be highly susceptible to orbital order be- 
tween the ions. Indeed, one of the aims of the present research is to study the influence 
of orbital ordering on optical transition strengths, and use this technique to study orbital 
ordering in YVO3. In this paper we present, amongst others, optical and magneto-optical 
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measurements of electronic excitations in YVO3 single crystals. The observed electronic 
spectra are discussed using ligand field theory considerations and compared to LDA+U cal- 
culations. In addition, magneto-optical Kerr effect experiments are used to elucidate the 
spin orientation. 

The crystal structure of yttrium orthovanadate YVO3 was initially reported to have an 
orthorhombic Phnm {D}^ in Schoenfiiess notation) symmetry at all temperatures^. Recent 
synchrotron measurements^*^ report the observation of an additional (401) refiection below 
200 K, inconsistent with the Phnm symmetry. The authors concluded^ that the crystal 
structure was most likely P2i/a, the highest symmetry subgroup of Phnm. However, the 
intensity of the Phnm forbidden refiection was four orders of magnitude weaker than to the 
allowed refiection^, which seemed to make the refinement of the X-ray data very difficult. 
Far infrared (FIR) spectroscopy is a very sensitive tool to study crystal structure variations. 
Minor atomic displacements due to a symmetry breaking, which may be difficult to elicit 
from x-ray analysis, manifest itself through the appearance of new phonon lines in the 
FIR spectrum. Recently, Ulrich et. al^ proposed that YVO3 can be subjected to a new 
sort of Peierls instability, where the new ground state involves a dimerization along the c 
axis. From the crystallographic and optical points of view this dimerization removes the 
inversion center, and IR and Raman active phonons become mixed. The a-axis modes stays 
nevertheless independent from the 6 or c axes, if the dimerization affects only the c axis. 
We use the FIR vibrational spectroscopy to follow temperature induced variations of the 
phonon spectrum and to refine the crystal structure of YVO3. 

II. EXPERIMENTAL DETAILS 

Polycrystalline powder of YVO3 was prepared from YVO4 powder as starting material. 
The YVO4 powder was prepared by a high temperature solid-state reaction from appropriate 
mixtures of pre-dried Y2O3 (99.998 %) and V2O5 (99.995 %, metal basis). The oxygen in 
YVO4 was reduced by annealing the powder in a fiow of pure H2 at 1000°C. The growth 
process was carried out by the fioating zone technique using a four-mirror furnace in a fiow 
of Ar gas^. The single-crystalline boule of about 6 mm in diameter and 60 - 70 mm in length 
was obtained from the growth. The crystallinity of the boule was checked by Laue X-ray 
diffraction. The elemental composition of the crystal was checked by the electron probe 
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FIG. 1: Pbnm orthorhombic crystal structure of YVO3. 



micro analyzer (EPMA) and the results showed that the molar ratios were given by Y:V:0 
= 1.00:1.00:3.02. A separate check of the composition using the chemical analysis method 
showed that the cation ratio of Y over V was 1.00 ± 0.01 and that the oxygen stoichiometry 
was 3.03 ± 0.02. Both results were in good agreement. 

The dielectric function for near-infrared, visible, and ultra-violent spectral ranges was 
obtained directly by using the ellipsometry technique. Because ellipsometric measurements 
are very sensitive to contamination of the surface, we used an ultra-high vacuum cryostat 
with a residual pressure less than 10~^ mbar. The cryostat was specially designed to minimize 
the movements of the sample due to the thermal expansion of the cold finger. The spectra 
were measured with a VASE32 ellipsometer from J. A. WooUam Co., Inc., covering a spectral 
range from 6000 cm~^ through 36000 cm~^ (0.75-4.5 eV). 

Far-infrared and mid-infrared reflectivity spectra were measured by Fourier Transform 
Infrared reflection spectroscopy using a Bruker IFS-113v. The intensity of the reflected light 
from the sample was referenced to that from a gold mirror for every measured optical surface 
at each polarization. The sample was mounted on the cold finger of a helium flow cryostat. 
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FIG. 2: The room temperature reflectivity and corresponding optical conductivity for the three 
crystallographic directions. Dashed lines indicate position of the phonon modes. There are 9, 9, 
and 7 modes for the a, b, and c axes, respectively, in agreement with the Pbnm symmetry. 

whose temperature was varied from 300 K down to 4.2 K. Kramers-Kronig transformation 
was used to calculate the dielectric function and optical conductivity from the reflectivity 
data in the infrared region. For these calculations we extrapolated the reflectivity below 50 
cm~^ with a pure ionic insulator response. For the high frequency extrapolation we used 
the reflectivity calculated from the ellipsometric data. 

The magneto-optical properties were measured using a home-built Kerr spectrometer. 
We used the electro-optical modulation technique, similar to [10], to obtain simultaneously 
the Kerr rotation and ellipticity induced by the sample. The combination of Xenon arc 
lamp, windows. Si-detector, and other optical components limited our spectral range to the 
band from 350 nm through 800 nm (1.55-3.5 eV). 
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FIG. 3: Temperature dependence of the FIR reflectivity for three polarizations. 
III. STRUCTURAL PHASE TRANSITIONS 

The orthorhombic Pbnm crystal structure of YVO3 is shown in Fig. ^ It can be con- 
sidered as a distorted perovskite structure. When compared to the perovskite unit cell, the 
orthorhombic a and b axes lay along the diagonals of the perovskite unit cell and the cor- 
responding translational vectors are approximately \/2 times longer. The c-axis dimension 
of the orthorhombic unit cell is twice as large. The YVO3 unit cell contains thus four for- 
mula units. Vanadium ions have an octahedral surrounding. The octahedrons are distorted, 
rotated and tilted with respect to each other, as illustrated in Fig. ^ 

In the present work we use optically active phonons to track changes in the crystal sym- 
metry. As yttrium orthovanadate is an insulating compound, its optical response in the 
far-infrared (FIR) region is determined purely by lattice vibrations. In Fig. |21 the room 
temperature reflectivity for all three polarizations is shown together with the optical con- 
ductivity. As can be seen, there are 9, 9, and 7 optically active modes respectively for the a, 
b, and c axes. The positions of the TO phonon frequencies are shown by vertical lines and 
the frequencies are given in Tables HI and ITTl Some resonances are very weak. For example, 
the mode of highest frequency for the c axis is hardly visible in the conductivity and appears 
as a shoulder on the reflectivity curve. It gains intensity as the temperature is lowered, and 
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FIG. 4: Left panels: "Leakage" of the b-axis i?3„ mode into the c-axis response. Right panels: 
Appearance of a new previously silent mode in the b and c responses. Insert: Temperature 
dependence of the optical strength of the new mode. 

reaches its maximum in the low temperature orthorhombic phase. Shown in Fig. 01 is the 
temperature dependence of the reflectivity for all three polarizations. The number of the 
phonons stays the same between room temperature and 200 K. Below 200 K a number 
of additional phonons gradually emerges. We illustrate this by picturing some of the new 
modes for the b and c axes in Fig. |3] As can be seen from the figure, the intensity of the 
new modes evolves continuously and monotonically as the temperature decreases, reaching 
its maximum at 75 K. The optical strengths remain however one or two orders of magnitude 
weaker than that of the original main modes. Below 75 K the new phonons discontinuously 
disappear. As an example, the temperature dependence of the oscillator strength of the new 
522 cm~^ mode is shown in Fig. EJ 

To analyze the observed phonon spectra, we apply a group theory analysis. The primitive 
cell of YVO3 contains four formula units, yielding 60 phonon modes in total. In the Pbnm 
lattice the four yttrium atoms and four vanadium atoms occupy equivalent sites 4c and 46 
in Wyckoff notation with the site symmetry C^^ and Cj, respectively. The 12 oxygen atoms 
are distributed between two inequivalent sites 4c and 8d with the site symmetry C^^ and 
Ci, respectively. Using the general procedure," we find that the phonon normal modes 
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FIG. 5: Transformations of the phonon irreducible representations due to the descent in lattice 
symmetry. 
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The total representation F can be subdivided into three acoustic vibrational modes, Biu + 
B2U + -Bsu, 24 Raman active modes, TA^ + 5i?ig + 7i?2g + 5i?3g, eight silent A„ modes, and 
25 infrared active modes. The infrared vibrations are split into 9-Bi„(£'| |a), 9-B3u(-E'| 
and 7B2u{,E\\c) modes with the electric dipole moment aligned along the a, 6, and c axes, 
respectively. 

A similar group theory analysis is applied to the lower symmetry structures, P2i/a, 
PI, and Phil, which are sub-groups of the Pbnm symmetry. Table ITTTl illustrates which 
symmetry elements are inherited in different symmetries. We also show a relation between 
the widely used notation Pbnm, adopted also in this paper, and the standard notation 
Pnma. The monoclinic P2i/a symmetry is postulated from the X-ray measurements^.. The 
triclinic PI symmetry, second highest possible symmetry after P2i/a, was proposed by us for 
the case, when the center of inversion is conserved^^. Both x-ray'*' and optical measurements — 
show that the symmetry can not be higher than monoclinic. Recent neutron experiments 
indicated the possibility of a dimerized state^i. Dimerization displaces the middle VO2 plane, 
see Fig. ^ away from its central position. It removes the inversion center and screw axis 
from the symmetry operations of P2i/a. The remaining symmetry group includes only a 



glide plane b and is referred as Pbll. 

As the X-ray measurements do not reveal any multiplication of the unit cell,^ the number 
of the atoms in the unit cell is not changed. For monoclinic P2i/a {€2}^ symmetry group one 
expects 16i?„ IR active and 2 acoustic modes in the he plane and 11 A,^ IR and one acoustic 
modes along the a axis. In case of triclinic PI {C}) symmetry, there should be SSA^ IR 
modes, which can have components along all the three axes. Monoclinic Phil {CD symmetry 
splits phonons into two irreducible representations 30A'(£'||6c) + 30^"(i?| |a), where Raman 
and IR modes are mixed. The way how the phonon modes' irreducible representations 
transform into each other due to the descent in symmetry is depicted in Fig. El Both 
monoclinic Cg^^ and triclinic C} crystal lattices have an inversion center. As a result the IR 
and Raman modes do not mix with each other for these symmetries. In point group C2h-i 
because the h and c axes are not orthogonal any more, the corresponding phonons mix and 
form the representation, 7i?2M + Q-Ba^ — IQBu- The mode A^, silent in D}^ symmetry, 
becomes optically active and together with the Biu mode gives the Au representation, 8A„ + 
9-Biu 11 Au- A possible twinning in the h — c plane cannot be a reason for leaking the 
6 — c phonons to the a axis and vice versa. Alternatively, in the lowest C} symmetry all the 
phonons mix together, as seen in Fig. |S1 In the Phil symmetry, Raman Ag modes become 
visible in the he plane and together with the IR By, phonons form the A' representation, 
l^Ag + 18i?u SOA'. Similar, Raman Bg modes can be observed in the a axis response 
together with IR Au modes, 12i?g + 18Au 30A". An important difference between PI and 
Phil symmetries is that in PI the same modes can be observed in all three polarizations, 
while in Phil a eigenmodes on one side and h and c eigenmodes on the other side remain 
separate. 

At room temperature we observed 9, 9, and 7 phonons for the a, 6, and c axes, respectively 
(see Fig. The TO phonon frequencies tu, broadening F and optical strength S are 
presented in Tables HI and HTl Our observation confirms the orthorhombic Phnm symmetry of 
the crystal lattice at room temperature. At lower temperatures, as far as the crystal structure 
is concerned, we can conclude from Fig. |3]that the Phnm symmetry is transformed to a lower 
one below 200 K and it is restored again below 75 K. Evolution of the optical strength of 
the new phonons with temperature indicates that the transition at 200 K. Surprisingly it 
is restored again below 75 K. The temperature evolution of the optical strength of the 
new phonons indicates that the phase transition at 200 K has a second order nature and 
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that the one at 75 K has a first order nature. Shown in Table |n] are the phonon mode 
frequencies and strengths observed in the c axis spectrum at 75 K. The lines active also in 
the high temperature phase are in bold face and assigned to B2U, those appearing only in 
the intermediate temperature phase are labelled by A. Although much more new phonons 
are seen for all three axes, only the activated c-axis modes, having a counterpart in the 
other axes, are shown in the table. First, one can see that the original c-axis phonon lines 
give rise to phonon absorptions along the b axis. The reverse is also valid. The phonons 
from the b axis 'leak' in the intermediate phase to the c axis, which is shown in Fig. |3J This 
agrees with all three symmetries under consideration. However, there is a number of newly 
activated absorption peaks (199, 350, and 522 cm~^) in the b and c polarizations, which do 
not have parents in other polarizations. The fact that these new modes are exactly at the 
same frequency in both b and c polarizations - as shown in the right panels of Fig. exclude 
the possibility of the doubling of the unit cell as their origin. The appearance of new modes 
is forbidden in the P2i/a symmetry (see also Fig. EJ. This observation rules out the P2i/a 
symmetry for the intermediate phase. Because these modes are not seen in any polarization 
above 200 K or below 75 K, we can ascribe these phonons either to the previously silent mode 
Au activated in the PI symmetry or to the high temperature Raman modes activated in 
Pbl 1 symmetry. We did not observe any correlation of phonons between b and c on one side 
and a on the other side, apart from the modes already present both in a and c polarizations 
above 200 K, 333 and 420 cm~^ modes in Table IHI This could imply that the triclinic 
distortions are extremely small, in particular for the a axis. However, the Pbll symmetry 
gives a more natural explanation, as it keeps the a and be modes apart consistent with the 
observation. All in all, these facts draw us to the conclusion that the crystal structure of 
YVO3 in the intermediate phase most likely has a monoclinic Pbll symmetry. A further 
test for this conclusion could be a comparison to Raman measurements, when data becomes 
available. 

Finally, we remark on the relation between the electronic and crystal structure changes. 
From Fig.inione can see that the main modes evolve smoothly through the second order phase 
transition, where the appearance of the orbital order is reported^. However, modes change 
their frequency substantially at the first order phase transition changes some, where spins 
and orbitals switch their ordering. The frequency change is found to be particularly strong 
for the modes around 400 cm^^. From lattice dynamical calculations^^ for the isostructural 
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FIG. 6: Temperature dependence of the TO frequencies for the c-axis i?2« phonon modes. For 
comparison purposes, vertical scale in each square comprises 6% of the central frequency. 

compound LaMnOs we conclude that the 400 cm~^ modes mainly involve vanadium-oxygen 
bond vibrations. Therefore, it seems that the main structural changes at the first order 
transition involve a distortion of the VOg octahedra. 



IV. ELECTRONIC STRUCTURE 

The yttrium Y^"*" and vanadium V^"*" ions in YVO3 both carry the formal charge 3+. 
Yttrium has a complete shell, and it can be excluded from the consideration for the electronic 
transitions. Vanadium has two electrons in the outer shell. We can thus expect two types 
of optical transitions in the visible range: d — d transitions between vanadium ions with the 
characteristic energy of the Mott-Hubbard gap (Mott-Hubbard (MH) transitions), and the 
charge-transfer (CT) transitions from the oxygen 2p-band to the free states in the vanadium 
3(i-band. 

Shown in Figs. [7| and |H1 are the real and imaginary parts of the YVO3 dielectric function, 
£ = £:' + %e" , measured at various temperatures. The optical response at low temperatures 
is determined by three bands with their maxima located at 1.8, 2.4, and 3.3 eV, followed 
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for YVO3. The vertical lines indicate locations of the band maxima. 
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FIG. 8: The real and imaginary part of the dielectric function for the a- and c-axes of YVO3. 
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by a sharp upturn above 4 eV. The optical gap is found at 1.6 eV. It is quite remarkable 
to observe that the intensity of some bands, located in the region of a few electron-volts, 
changes by as much as 50% for temperature variations as low as fc^T < 26 meV. 

Recently Miyasaka et. al}^ reported reflectivity measurements on YVO3, where the op- 
tical conductivity, obtained through the Kramers-Kronig relations, revealed only two peaks. 
Also the temperature dependence and the reported anisotropy of the conductivity were dif- 
ferent. Since e" is approximately five times smaller than e' in the visible part of the spectrum, 
the reflectivity is mainly determined by the real part of e. The rf-band structure in e' and 
e" gives rise to the corresponding reflectivity variations of only few percent and the temper- 
ature induced changes of the reflectivity are around 1% or less. This together with the close 
vicinity of the strong charge transfer transitions place heavy demands on the accuracy of 
the measurements. A small error can substantially affect an interpretation of the electronic 
spectrum. We used the ellipsometric technique, which determined e' and e" directly at each 
particular photon energy and excluded the necessity of reference measurements. In addition 
we took special precautions to perform measurements in ultra-high vacuum in order to avoid 
condensation on the surface at low temperatures. We trust that these measures minimize a 
possible error. 

Let us first discuss the origin of the bands. Following the general classification scheme^^, 
YVO3 is a Mott-Hubbard insulator. The lowest excitations are dfd?- d^d}- transitions 
between different vanadium ions, i and j. The upturn above 4 eV is due to charge transfer 
excitations. This general picture follows from our LSDA+f/ calculations and agrees with 
previous room temperature measurement o^^ 1 . The fine structure of the MH band observed 
here has not been reported previously and will be discussed using ligand field theoretical 
considerations and LSDA+U calculations. 

Because the orthorhombic distortion of the oxygen octahedron is small, we can neglect 
any on-site d-d excitations. For the same reason we can, considering the energy diagram, 
restrict ourselves in the first approximation to the Oh symmetry on the vanadium site. The 
eg orbitals can be excluded, because due to the large crystal field splitting, lODg, the mixing 
between Cg and t2g orbitals is expected to be small. Thus, in the strong field approximation 
the two d electrons occupy t2g orbitals. 

If we take into account the Coulumb interaction, the t2g level is split into '^Ti^, ^Eg, 
^T2g, and ^Aig levels, as shown in Fig. IHl (a). The ^Eg and ^T2g levels are degenerate in 
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FIG. 9: Energy diagram for (a) the ground state of the vanadium ions and (b) excited states of 
two vanadium sites. 

the Oh symmetry. In the ground state the d electrons occupy the lowest energy S=l state, 
^Tig, in agreement with the Hund's rule. The MH excitations occur from the initial state 
with two vanadium in the d'^ state to the final state with the first atom in the d^ (V*^^) 
and the second in the d^ (V^'^) states. In the final state the wave function of the single 
electron follows obviously only one irreducable representation '^T2g, while the d^ state splits 
into four levels, '^A2g, "^Eg, ^Tig, and ^T2g, with "^Eg and ^Ti^ again being degenerate. The 
energy diagram of the final state is depicted in Fig.|ni(b). The three observed bands in Fig.EI 
correspond to transitions to these three energy levels. According to the measurements, these 
bands span the energy range of 1.5 eV from 1.8 eV to 3.3 eV. With the crystal field splitting 
of lODq ^ 2 eV, our assumption of a small admixture of Cg orbitals is justified. 

The picture obtained with the ligand field theory is in good agreement with our LSDA+t/ 
calculations. The LSDA+[/ band structure calculations has been performed in the LMTO 
calculation schem o^^i^^ . Vanadium 3d partial density of states were obtained in the LDA+U 
calculation with U = 3.4 eV and J = 0.85 eV for the low temperature crystal and magnetic 
structure of YV03. The results are consistent with previous calculations of the electronic 
structure^'^. As can be seen from the electron density of states in Fig. the highest 
occupied states are from a narrow d band, with both electrons in a spin-up state. The AFM 
ordered neighboring vanadium ion should be accordingly in a spin-down state. The lowest 
unoccupied band corresponds to the addition of an electron with spin up. This agrees with 
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FIG. 10: Vanadium 3d partial density of states for spin-up and spin-down electrons, obtained in 
local spin-density approximation with the Hubbard potential U taken into account. — 1 and 
-|- 1 correspond to the electron removal and electron addition states, respectively. 

the lowest excited '^A2g state, which indeed puts electrons in the high-spin state S = 3/2. 
Two consequent excitation bands are formed by transfer of spin-down electrons, again in 
agreement with the low S = 1/2 spin states in the "^Eg, "^Tig, and '^T2g representations found 
from the ligand field theory. Higher electron-addition bands with equal spin- up/spin-down 
populations have an Cg origin. From the energy distribution of the bands in Fig. one 
can see that YVO3 is indeed a Mott-Hubbard insulatop^^. The lowest excitations are d — d 
transitions, d'^d'^ — > d^d^ . The charge-transfer transitions, rf^ — > d^L occur at higher energies 
than the d — d transitions. 

The temperature dependence of the optical response of YVO3 is shown in Figs. [7| and 
IHl There is only small T-dependence between 70 K and 6 K, as can be seen in Fig. |Hl 
and the corresponding data are further omitted from the consideration. To clarify the 
temperature dependence of the MH bands. Fig. ^2 shows a differential dielectric function 
for E\\h^ where the 250 K data are subtracted from the curves in Fig. [7| In order to quantify 
temperature variations of the optical strength of the bands, we fitted simultaneously the 
real and imaginary parts. The differential dielectric function turned out to be well fitted 
with three gaussians in the imaginary part and with its Hilbert transform, Dawson function. 
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FIG. 11: The dielectric functions from Fig. [71 the real and imaginary parts, with subtracted 250 
K data. The fit curve shows a simultaneous fit to the real and imaginary part for the 80 K data. 

in the real part. The fitted curves for 80 K are shown along with the experimental data 
in Fig. ^2 The full dielectric function shown in Fig. [7| was also fitted in the same way, 
but an additional fourth gaussian was necessary to account for the CT band. In this way 
the position of the bands was confirmed to be temperature independent. However, it was 
not possible to determine accurately the absolute optical strength of the MH bands due to 
an additional non-gaussian contribution. It was also impossible to determine the position 
of the 3.3 eV band at high temperatures because of its vanishing presence in the spectra. 
The band positions and optical strengths obtained from the fits are plotted in Fig. ^| The 
dielectric function reveals minimal variations between room temperature and the second 
order structural phase transition at 200 K. Below 200 K a strong growth of the 3.3 eV 
band is observed, which is accompanied by the decrease of the 1.8 and 2.4 eV bands. The 
total change of the optical strength of these MH bands remains virtually zero down to 
the Neel temperature, Tn = 116 K. In the magnetically ordered state the 3.3 eV band 
acquires additional strength. The additional spectral weight transfer should come from 
higher energies, most probably from the Cg states, as the decrease of the low energy bands 
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FIG. 12: Temperature dependence of the bands position (left) and their optical strength together 
with the total spectral weight of the three bands (right). 

cannot account for it. Indeed, the lowest energy state in the t^gCg configuration is obtained 
by adding an eg electron to the ground state ^Tig. Four irreducible representations, "^Tig, 
^T2g, ^Tig, and ^T2g can be formed from this configuration, where the latter representation 
descends also from the tf^ configuration. Antiferromagnetic order favors transitions to the 
low-spin states and thus gives additional optical strength to the '^T2g state. 

At the first order phase transition, the optical strength of the 2.4 and 3.3 eV bands 
experiences jumps in opposite directions. Additionally, as can be seen in Fig. ^2 (a), a new 
shoulder emerges at around 2.7 eV. This appears in Fig. Has an extended flat region, which 
becomes visible probably due to the narrowing of the 3.3 eV band. However, its optical 
strength is very small. 

The optical spectrum begins to change below the second order phase transition at 200 K. 
These changes can not be directly related to the lower symmetry as the 3.3 eV band does 
not disappear in the low temperature orthorhombic phase. The reason should rather come 
from some electronic correlations and associated distortions of the VOe octahedrons. Ren et 
ali^ suggested that the spin order in YVO3 is accompanied by an orbital order. The orbital 
order is reported to appear below the 200 K phase transition^. In the proposed picture for 
the orbital ordering, the xy-orbitals are always occupied by one of the two d electrons, while 
the occupation of the xz and yz orbitals alternates in an AFM order. The orbital order was 
found to be different above and below 77 K in agreement with the different magnetic order. 
According to Ren et al^, there is a C-type magnetic and G-type orbital order above 77 K, 
and alternatively G-type magnetic and C-type orbital order below 77 K. 
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We now turn to the discussion of how orbital and magnetic ordering influences the op- 
tical conductivity. We consider a model, where the optical conductivity is proportional to 
the electron transition probability between two neighboring vanadium sites. There are in 
total nine independent eigenfunctions in the ^Tig representation. At room temperature the 
electron transition probability can be calculated by averaging over arbitrary orientations of 
spin and orbits on two ions, i.e. by averaging over all linear combinations of these 9 wave 
functions. Magnetic or orbital order limits the choice of the initial wave functions. At lower 
temperatures we can expect that the electron transition from the S=l state on the first ion 
to the S=3/2 state on the other ion becomes less probable due to the antiferromagnetic 
order. Thus, the intensity of the lowest band should decrease. As can be seen from Fig. [7| 
and Fig. ^2 the optical strength of the lowest band indeed decreases together with tem- 
perature. By using a more rigorous analysis of the transition probabilities, we can evaluate 
if a particular order is consistent with the observed temperature dependence for e{uj). In 
this analysis we neglect the tilting of octahedrons, and assume that the transition integrals 
between rf-orbitals on the neighbor vanadium ions are the same for all directions and for all 
d orbitals. The conductivity is assumed to be determined by the probabilities to create a 
pair consisting of an occupied orbital on one site and an empty orbital on the neighboring 
site, with the constraint that the transition is allowed. The intersite d — d transitions are 
mediated by the p orbitals of the oxygen ions between the vanadium sites. At room temper- 
ature all ^Tig states are occupied randomly. At lower temperatures we used the magnetic 
and orbital order, suggested above^. If we assume the transition integral from completely 
occupied to completely empty orbital as unity, then the transition probabilities at room 
temperature are found to be 0.49, 0.84, and 0.44 for the ^^29, "^Eg +^ Tig, and ^T2g final 
states, respectively. We also calculated the transition probabilities for AFM spin order and 
for AFM spin order accompanied by the aforementioned orbital order (see Table HV|) . Both 
simple AFM- SO and AFM-SO+00 lead to a decrease of the lowest band and an increase 
of the intensity of the highest band, consistent with the observations. 

However, none of the considered orderings gives the correct sign for the intensity change 
for the middle band, i.e. a decrease of strength with decreasing temperature. Moreover, we 
found that no combination of completely polarized states {Sz = 1) can give the desired in- 
tensity decrease for the "^Eg +^ Tig band. The same conclusions hold if the on-site symmetry 
is lowered to tetragonal. If some addition of the Sz = initial wave function is allowed, as 
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FIG. 13: (a) Temperature dependence of the Kerr rotation measured in the polar geometry with 
H\ \a. (b) The Kerr rotation Qk together with the Kerr elhpticity ex in the low temperature phase. 
(Insert) Enlarged view of the resonances at 1.8 and 2.2 eV. 

suggested by neutron measurement s4, a whole number of possible orbitally ordered struc- 
tures becomes allowed. Although different orbital orderings can indeed explain the observed 
intensity changes, it is not possible to give a preference to a particular order—. 



V. MAGNETO-OPTICAL KERR EFFECT 

The energy diagram shown in Fig. El was derived assuming a local cubic symmetry on 
the vanadium site. In reality, one axis of the octahedron is different from the other two. 
This tetragonal type distortion additionally splits the multiplet structure, which partially 
lifts the degeneracy of the ground and excited states. On top of that, a spin-orbital coupling 
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provides segregation of the states by the quantum number Mj. The latter gives rise to 
the magneto-optical Kerr effect (MOKE). Temperature dependence of the Kerr rotation 
spectrum Qk{E) is shown in Fig. The curves are shifted for clarity. The measurements 
were performed in the polar geometry at nearly normal angle of incidence, with magnetic 
field of 0.1-0.2 T parallel to the a axis, H\\a. In the other field orientations, the magnetic 
moment was too small for the Kerr rotation to be detected. No rotation of the polarization 
plane was observed above the Neel temperature within our experimental error. Below the 
Neel temperature but above the first order phase transition, only a weak structure of the 
order of Qk ~ 0.01° can be seen. Below 75 K a remarkable peak as large as = 0.1° is 
observed. In total, there are three resonance structures. Enlarged in the insert of Fig. 
are two surprisingly narrow dispersion type resonances at 1.8 and 2.2 eV that have positions 
approximately corresponding to the onset of the two first MH bands. Their intensity, though 
rather small ~ 0.01°, is well reproduced on different samples. The third broad peak appears 
to have maximum at 3.5 eV, which is also confirmed by the zero crossing of the Kerr 
ellipticity 6k{E). This band lies below the charge-transfer onset and corresponds to the 
d — d transitions. The maximum rotation in this band amounts to ~ 0.1°. Such a large 
magneto-optical effect seems to be quite surprising, because from the electrodynamics point 
of view the time-reversal symmetry, although broken on the microscopic scale, is usually 
restored on the macroscopic scale in antiferromagnetic compounds. If this indeed were the 
case in YVO3, the Kerr rotation would mainly result from the small 0.2° spin canting and 
would be as minuscule as the net magnetic moment is. For comparison, the saturation 
magnetic moment in Ni is 0.61/2^ per atom, i.e. approximately 60 times larger that the net 
moment of 0.01 Bohr magneton per vanadium atom in YVO3. Still the Kerr rotation in Ni 
is quite comparable to those of the low temperature phase of YVO3 and lies between 0.1° 
and 0.2° in the visible range. 

Antiferromagnetic compounds can have large magneto-optical effects, as was shown 
recent ly2S>2i, if the local spin and orbital magnetic moments are not collinear. In partic- 
ular, an orbital magnetic moment does not vanish in transition-metal perovskites with the 
Phnm structure due to the tilting of the octahedra^i. We can find non-vanishing ferromag- 
netic components of the net orbital magnetic moment from symmetry considerations, 
following the guideline and notations from Ref. |3 There are four inequivalent vanadium 
sites i = 1,2,3,4 in the YVO3 unit cell. These sites can be generated by four symmetry 



21 



transformations, unity operation E, and three screw axes, S2x with the shift (|,0,0) and 
coordinates x,i,0; S2y{0,^,0) ^,y,^; and 5*22(0,0,1) 0,0, z. Here x,y, and z are parallel to the 
crystallographic axes a, 6, and c, respectively. The local relation between the orbital mag- 
netic moment and spin magnetization direction e = (cos 6 sin 0, sin 6 sin 0, cos 6) on each 
vanadium site i can generally be written as = £^^e/3, where the matrix determines 
the orbital magnetic moment induced by the spin-orbit interaction. We apply the mentioned 
symmetry operations to the matrix to obtain averaged matrices Cais = J2t=i ^^ap 
the C- and G-type spin ordering, respectively: 



^ ^ 



c 



\ 








c 



G 



CI 




CI 



The corresponding net orbital magnetic moments are 

C : Ml = 4:{Cly sine sin (f),Cly sine cos (f),0), 
G: = 4{Cl^cose,0, Cl^ sine cos (p). 

The orbital magnetic moment has the same symmetry as the gyration vector g2^, which 
determines the antisymmetric part of the dielectric tensor, = ieap-yg^^, which is respon- 
sible for mageto-optical effects^^™. Here e^/j^ is the unit antisymmetric tensor. In the polar 
measurement geometry discussed here, components of the vector Qa and correspondingly 
components of define the orientation of magnetic field H , at which the observation of 
Kerr effect is possible. 

For the G-type spin order, the structure of suggests that the polar Kerr effect can 
exist for H\\a, if spins are parallel to the c axis, and for H\\c, if spins are aligned along the 
a axis. If the spins are along the h axis, and the corresponding Kerr effect are zero. 
This does not exclude the possibility for the second, weaker, type of Kerr effect, which is 
directly related to the spin canting. It should be emphasized here that the Kerr rotation due 
to the orbital magnetic moments can exist even if spins are exactly antiferromagnetically 
aligned^i. Our observation of the substantial magneto-optical effect only for H\\a shows 
that spins in the low temperature phase of YVO3 are aligned along the c axis. This agrees 
with the neutron scattering measurements^. The weak structure at 1.8 and 2.2 eV is most 
probably the spin canting effect. 
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The C-type spin order as seen from the M£ angular dependence allows a magneto-optical 
activity, only if spins are aligned in the ab plane. The small magnitude of the experimentally 
observed MOKE above 80 K compared to the low temperature phase suggests that the spins 
are aligned along the c axis. In this case, M£ is zero for any orientation of H, and Kerr 
effect substantially decreases, in agreement with the our experimental observations. This 
result contradicts the neutron data, where a substantial spin component along the b axis is 
observed^. This, however, can reflect another important property of YVO3. As it was shown 



in Ref. 



271 the G-type antiferromagnetic orbital order in the form proposed by Ren et. al^ 



suppresses the orbital magnetic moment. This can be the reason, why no substantial Kerr 
rotation was observed for the C-type spin order. This conclusion requires that the on-site 
magnetic moment in the G-type orbital order is reduced from 2^b to not due to 

the orbital magnetic moment, but due to a mixture oi = 1 and 5*2 = ^Tig states, in 
agreement with the conclusion obtained in the previous section. 



VI. CONCLUSIONS 



The analysis of the YVO3 phonon spectra revealed three different phases. The orthorhom- 
bic Pbnm structure exists down to 200 K. At 200 K there is a second order phase transition 
to a lower symmetry structure. Orthorhombic Pbnm crystal symmetry is restored through 
a first order phase transition below 74 K. The most probable structure in the intermediate 
phase is the monoclinic P611, where the VOg octahedra form a dimerized chain along the c 
axis. The triclinic PI symmetry cannot however be completely excluded. 

YVO3 is a Mott-Hubbard insulator with the optical gap of 1.6 eV. There are three optical 
bands with energies 1.8, 2.4, and 3.3 eV. We identify these bands with the inter-vanadium 
transitions from the ^Tig ground state to the ^^29, "^Eg Tig, and ^T2g states. These 
d — d excitations are followed by the charge-transfer transitions at 4 eV. The band gap 
structure is in agreement with the LSDA+f/ calculations. The d bands revealed strong tem- 
perature dependence. Using a model, which includes only the nearest neighbor transitions, 
we found that the observed temperature dependence disagrees with the simple antiferro- 
magnetic spin ordering. This indicates that an orbital order is responsible for the band 
temperature variations. However, the commonly considered \xy,xz\/\xy,yz\ orbital order 
seems also inconsistent with our data, if we assume a fully polarized on-site state. 
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We observed two types of magneto-optical effect. Small net magnetic moment due to 
the uncompensated spin moments induces the weak Kerr rotation of the order of ~ 0.01°. 
The effect exists in all phases below the magnetic ordering temperature and for all active 
optical transitions. A strong Kerr effect appears as a consequence of the low crystallographic 
symmetry and ferromagnetic ordering of orbital magnetic moments. The Kerr rotation of 
the order of ~ 0.1° was found in the low temperature phase for the highest d — d transition 
band. It follows from the symmetry considerations that the spin magnetic moment lays 
along the c axis in both magnetic structures. However, this conclusion for the C-type AFM 
spin arrangement can be affected by the accompanying orbital order. 
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TABLE II: c-Axis phonons at 75 K. modes originate from the room temperature phase. New 
vibrational modes, appearing below 200 K, are designated as A 



Mode 


Frequency 


Sa 




So 




180.2 


— 


— 


7 


A 


199 




0.037 


0.008 


B2u 


205.6 




0.006 


0.78 


B2u 


333.0 


1). If) 


0.008 


2.9 


A 


350 




0.67 


0.003 


B2u 


359.8 




0.02 


0.37 


B2u 


420.5 


1.6 


0.076 


1.6 


A 


522 




0.024 


0.018 


B2u 


564.9 




0.014 


1.2 


B2u 


699.2 






0.002 



TABLE I: TO phonon modes for the a and b axes at room temperature. 





Biu (a) 






B3u (b) 




u, cm ^ 


7, cm~^ 


Sa 
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7, cm~^ 


Sb 


162.2 


3.2 


0.08 


144.3 


15.5 


0.73 


211.6 


8.5 


5.4 


237.4 


5.5 


2.0 


328.5 


9.9 


1.6 


280.1 


8.4 


0.64 


359.2 


15.0 


0.56 


318.7 


7.8 


0.35 


420.8 


10.4 


1.4 


374.5 


6.3 


0.22 


446.8 


17.3 


2.6 


390.7 


16.5 


4.21 


494.5 


10.0 


0.005 


471.7 


13.3 


1.2 


517.5 


15.5 


0.06 


543.7 


9.8 


0.86 


572.4 


22.6 


0.91 


559.8 


17.5 


0.21 
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TABLE III: Possible symmetries and symmetry operations in the intermediate phase {75K < T < 
200K). 
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Pii run 


Phniii 


.P2i/« 


PI 


Pbll 




r)i6 


r)l6 




CI 
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2(0,0,i)i,0,z 


2(i,0,0)x,i,0 


2i 


- 
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2(0,i,0)0,y,0 


2(0, 0,i) 0,0,2; 








(4) 


2(i,0,0)x,i,i 


2(0,i,0)i,2/,i 








(5) 


I 


z 








(6) 






b 




b 


(7) 


mx, \,z 


mx,y, I 








(8) 


n(0.1.i)i.//.^ 


n(i.0.i),r.i.. 









TABLE IV: Calculated transition probabilities at room temperature, in the AFM spin ordered 
(SO) state and in the AFM spin and orbital ordered (SO-00) state. 

Order Tig 

Room temperature 0.49 0.84 0.44 

AFM SO 0.30 1.48 0.89 

AFM SO-00 0.33 1.67 1.0 
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